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TEOREMAS DE CONVERGÊNCIA NOS MÉTODOS VARIACIONAIS (*) 


RESUMO 


Apresenta-se uma teoria dos métodos variacionais 
a qual pressupõe a existência de um princípio de mi- 
nimo. 

A solução exacta do problema variacional supõe-se 
pertencer a um espaço de Banach, X, onde cada ele- 
mento minimiza um determinado funcional numa certa 
classe de elementos. 

Considera-se um operador que faz corresponder a 
cada elemento de X um elemento de um dos subespa- 
cos de X, o qual elemento se denomina aproximação 
do primeiro. 

Enunciam-se e demonstram-se teoremas que esta- 
belecem limites superiores para a distância de cada 
elemento X e a sua aproximação. 

A teoria pode ser utilizada em situações mais ge- 
rais que a correspondente ao método de Ritz. Tem, 
muito especialmente, o interesse de permitir enqua- 
drar os problemas de convergência do método dos ele- 
mentos finitos. 


1 — INTRODUÇÃO 


Os problemas de Física Matemática reduzem-se 
geralmente ao da resolução de equações diferenciais às 
quais está muitas vezes associado um princípio varia- 
cional [1] que afirma que a solução relativa a certas 
condições de fronteira minimiza um certo funcional 
num espaço de funções. 

Os métodos variacionais são justamente aqueles 
que substituem o problema da resolução da equação 
pelo da determinação do elemento que minimiza o 
funcional. 

O presente trabalho considera o problema já na 
fase em que ele se apresenta como problema variacio- 
nal. Estabelecem-se os teoremas gerais de aproxima- 


(*) Comunicação apresentada nas 1.ºº Jornadas Matemáticas 
através do projecto TLE4 do Instituto de Alta Cultura. 
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C.D.U. 519.34 
EDUARDO ROMANO DE ARANTES E OLIVEIRA 
Professor Catedrático do Instituto Superior Técnico, Lisboa. 
SUMARY 


A theory of variational methods is presented which 
supposes a minimum principle. 

The exact solution of the variational problem is 
supposed to belong to a Banach space, X, where each 
element minimizes a certain functional on a certain class 
of elements. 

An operator is considered which associates to 
each element of X an element belonging to a subspace 
of X. This element is called the aproximation of the 
first. 

Theorems are stated and demonstrated which esta 
blish upper bounds for the distance between each ele- 
ment in X and its approximation. 

The theory can be used in more general situations 
than the one corresponding to Ritz's method. It provides 
namely an ideal frame for analyzing convergence in the 
finite element technique. 


ção e convergência relacionados com a redução do pro- 
blema variacional inicial a outros problemas também 
variacionais mas de resolução eventualmente mais fácil. 


2 — Subconjuntos isocondicionados e isominimizantes 

Seja H um espaço de Hilbert real [2] e seja F 
uma família de funcionais contínuos definidos sobre H. 

Supõe-se que cada funcional da família F admite 
um único elemento minimizante global num subespaço 
He € H, dito um subespaço condicionado de H, e numa 
classe de subconjuntos de H, chamados subconjuntos 
isocondicionados, cada um dos quais se obtém somando 
um elemento constante aos elementos de H.. 


Luso-Espanholas, Lisboa, Abril de 1972. Trabalho subsidiado 
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isto significa que cada subconjunto isocondicionado 
Cc H, contém, para cada funcional F, um elemento s 
tal que 


Fis) < F(c) (2.1) 


sendo c um elemento arbitrário de C distinto de s. 


Cada subconjunto isocondicionado C e o subes- 
paço H. são isométricos, o que assegura que cada 
elemento de C seja um ponto de acumulação [2] de C. 
É além disso evidente que o subespaço condicionado é 
um dos subconjuntos isocondicionados e que a união 
de todos estes é o espaço total H. 


Considere-se agora, para cada funcional F da famí- 
lia F, o conjunto de todos os elementos minimizantes 
contidos em todos os subconjuntos isocondicionados. 
Designa-se este novo subconjunto por subconjunto iso- 
minimizante correspondente a F. Supõe-se que a família 
F é tal que a união de todos os subconjuntos isomini- 
mizantes, correspondentes a todos os funcionais da fa- 
mília F, coincide com o espaço todo, H. Cada ponto 
de H representa pois a intersecção de um subconjunto 
isocondicionado com um subconjunto isominimizante. 


Dois elementos pertencentes ao mesmo subcon- 
junto isominimizante dir-se-ão elementos isominimi- 
zantes. Dois elementos pertencentes ao mesmo subcon- 
junto isocondicionado dir-se-ão elementos isocondi- 
cionados. 


Os subconjuntos isocondicionados e isominimi- 
zantes supõem-se ortogonais [3], Por isso, a intersec- 
ção de cada subconjunto isocondicionado com cada sub- 
conjunto isominimizante, que se supõe sempre não 
vazia, não pode conter mais de um elemento. 


Considere-se agora um operador limitado e conti 
nuo, T, com domínio H e contradomínio H”, sendo H' 
um subespaço de H. A cada elemento de H corresponde 
um único elemento de H' mas a cada elemento de H' 
corresponderá em geral mais de um elemento de H. 
Supõe-se que a imagem de qualquer elemento de H per 
tencente a H' coincide consigo mesmo, isto é, que 


T(e)=e se ec H' (2.2) 

Supõe-se ainda que podem definir-se em H' sub- 
conjuntos isocondicionados e isominimizantes, relativos 
à mesma família de funcionais F, subconjuntos esses 
que gozam das mesmas propriedades que os subcon- 
juntos homólogos em H. 


As imagens dos elementos de cada subconjunto 
isocondicionado C de H supõem-se contidas num sub- 
conjunto isocondicionado de H', o qual se diz corres- 
ponder a C. Isto significa que as imagens de quaisquer 
dois elementos isocondicionados, e e f, de H são dois 
elementos isocondicionados, e' e f”, de H (fig. 1), em- 
bora não sejam necessáriamente elementos isocondicio- 
nados de H. Se no entanto o forem, quaisquer que sejam 
e e f, o operador T diz conformante. (*º) 


e tn 
Fig. 1 


É evidente que, se um dado elemento pertencer a 
H', então, em virtude de (2.2), coincidirá com a sua 
própria imagem, Estará pois contido simultâneamente 
em Ce em C (fig. 2). 


Fig. 2 


Um certo subconjunto isominimizante M de H' 
diz-se corresponder a um certo subconjunto isominimi- 
zante M de H se ambos corresponderem ao mesmo 
funcional, F. 

Pode pois considerar-se um segundo operador, A, 
com domínio H e contradominio H, o qual faz corres- 
ponder à intersecção de cada par de subconjuntos iso- 
condicionados e isominimizantes de H a intersecção 
dos subconjuntos correspondentes em H'. Tendo em 
vista o que ficou estabelecido atrás, o operador A é 
também tal que 


Alej=ze se ee H' (2.3) 


Exactamente como acontece com os subconjuntos 
isocondicionados, a intersecção de qualquer subcon 
junto isominimizante de H com o correspondente sub- 
conjunto isominimizante de H' é um elemento de H' 
que coincide com a própria imagem (fig. 2). 

Sobre cada elemento de H' cruzam-se pois um 
subconjunto isocondicionado e um subconjunto isomi- 
nimizante de H, e os subconjuntos correspondentes 
de H. 


(*) Esta designação é inspirada pela linguagem utilizada na teoria dos elementos finitos [4]. 
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À A — imagem de e dá-se o nome de elemento 
aproximante ou aproximação de e em H.. 


3 — Os teoremas de aproximação 


Os teoremas de aproximação apresentam condi- 
ções suficientes para que a aproximação de um dado 
elemento de H caia dentro de uma certa vizinhança 
do elemento. 

Seja s o elemento minimizante único de F no sub- 
conjunto isocondicionado C de H. Seja s - o elemento 
minimizante único de F no subconjunto isocondicio- 
nado correspondente, C', de H' . s. é assim a apro- 
ximação de ss. 

Considere-se em primeiro lugar o seguinte lema: 


c 


Lema: «Dado um número positivo arbitrário 5, é 
possível determinar outro número positivo arbitrário 5' 
tal que, se c for um ponto contido em C e F(c) — 
— F(s)< 5º, a distância entre c e s será mais pequena 
que 5 ». 

Considere-se com efeito uma esfera Sy em C, 
com raio ) e centros, e seja Fy o valor minimo que 
F toma no conjunto C-S 3 , isto é na porção de C ex 
terior à esfera. Em virtude de (2.1), 


F3-F(s)> O (3.1) 
Assim, dado qualquer ponto C tal que 
F(c)< F3 (3.2) 


a distância d (s,c) entre s e c deve ser mais pequena 
que 5, o que significa que 9%" pode ser igual a 
F3x -Fd(s), isto é que 3 pode ser a diferença entre os 
valores mínimos tomados por F em C, fora e dentro 
da esfera Sa. 

Um segundo lema poderia ser imediatamente esta- 
belecido adaptando o anterior ao espaço Hº'. 


O primeiro teorema de aproximação é o seguinte: 


«Dado um número positivo arbitráriamente pequeno, 
5, é possível determinar um outro número positivo arti- 
trãriamente pequeno, :, tal que, se 


a) o operador T for conformante e C' estiver incluído 
em C, 
b) existir um ponto c' E C' cuja distância a s seja 
menor que :, 
então a distância d (s, s',) é menor que 5». 
A figura 3 traduz graficamente a situação. 


Fig. 3 
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Como C' CC e SEC, 


F()<F(s,) (3.3) 
Como cel 
F (s,)< F(c) (3.4) 
Em virtude da condição b), 
disc)< e: (3.5) 


Ora, como F é contínuo sobre H, é possível, dado 


um número positivo arbitrariamente pequeno, +: , deter- 
minar : tal que (3.5) implique 
|IF(c)— F(9)|< +! (3.6) 
Combinando (3.3) e (3.4), obtém-se 
F()< F(sS,)< F(e) (3.7) 
e, em virtude de (3.6), 
F(s)< Fis )<F()++ (3.8) 
o que implica 
FP (s)— F(s)<s (3.9) 


Como s, pertencea C e C CC é então possí- 
vel, em virtude do Lema, determinar :«, e portanto :, 
de tal modo que 


dis s,)<3 (3.10) 


Considere-se agora o segundo teorema de aproxi- 
mação. 


«Dado um número positivo arbitráriamente pequeno, 
3, é possível determinar dois outros números positivos e 
arbitáriamente pequenos : e :' tais que, se: 


a) existir um ponto c' E C' cuja distância a s seja 
menor ques, 
b) existir um ponto c& C cuja distância a s , seja 
menor que c', 
a distância entre s e s', é menor que 5 ». 
A situação é descrita gráficamente pela fig. 4. 


Fig. 4 
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Com efeito, como ce €C, 


F(s)< F(c) (3.11) 
e, como c E C”, 
Fis, )< F(c') (3.12) 
Em virtude das condições b) e c), 
d (s,c') (3.13) 
o (6,8, | =0*" (3.14) 


e, como F é um funcional contínuo sobre H, é possível 
escolher : e :* de tal modo que 


F(s)=F(c)+:” (3.15) 
Fl)=Fis,)+te” (3.16) 
3 (3.17) 


sendo : “ e : '* números positivos arbitrariamente pe- 
quenos. 

Eliminando Fis", ) e F(c') entre (3.12), (3.15) 
e (3.16), e usando (3.11), obtém-se 


Fig <Fio<FisrE rr (3.18) 


e portanto 
Flc)-—F(g)< e” 4 1º (3.19) 


Em virtude do Lema, é pois possivel determinar 
“e portanto : e: tais que 


em 


dics)< 53 —:s' (3.20) 


A desigualdade triangular permite então escrever: 


dis,s',)< dfc,s) d(c,s',) (3.21) 
e portanto, em virtude de (3.14) e (3.20), 
dis.s',)<s (3.22) 


Nas aplicações destes teoremas, convém que o 
elemento que em ambos os teoremas se supõe existir 
na vizinhança : do s seja a T — imagem de s, isto é, 
T (s). 

Quanto ao elemento c & C que, no segundo teo- 
rema, se supõe localizado na vizinhança :' des, deve 
dizer-se que, na prática, não é fácil mostrar a sua exis- 
tência. 

Pode ser conveniente utilizar como elemento c um 
elemento s, de C tal que T(s,)= Ss", - Outra pos- 
sibilidade consiste em tomar para elemento c a projec- 
ção p de s', sobre C (fig. 5), recorrendo ao seguinte 
teorema: 
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«Seja p a projecção de s', sobre C e p' a 
T — imagem de p em H'. Seja z a norma (suposta 
e! 
limitada) do operador A. Então, se d (p, p') ' 
+ 2 
a distância d (s',, p) é menor que <'». 


Mostra a fig. 5 que p' e s', estão ambos em 
C'. p e q são projecções de s', e p' sobre C. 

Pelo teorema de Pitágoras, válido para espaços 
de Hilbert; 


t ut 49 
[d(p.q)]2 [d(p'q)]2 = [dtp,p')]?2 E?) TI, 


(3.23) 


Fig. 5 


Portanto, 


(3.24) 


Mas M , [5 M'q intersectam-se sobre o elemento 
p' que pertence a H'. Portanto, M' é o subconjunto 
isominimizante de H' que corresponde a M a O mesmo 
pode ser dito de M , º M' visto que eles se inter- 
sectam sobre s',. 

Como, por outro lado, p' e s', estão contidos em 
C', que é o subconjunto isocondicionado de H' que 
corresponde a C, segue-se que p' e s', são as apro- 
ximações de q e p. Logo, 


p + 


d(p',s',)< a dí(p,q) (3.25) 
sendo z a norma do operador A, 
Ora 
d(s',.p) <. d(s',, p) + d(p,p' (3.26) 
e, portanto, em virtude de (3.24), 
d(s',.p) «+ (3.27) 
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4 — Convergência 


Consideram-se agora sucessões de operadores limi- 


tados e contínuos, [ La h e ( A, ) - À cada termo 
destas sucessões corresponderá um subespaço H'. C H 
que é o contradomínio dos operadores T, e A,. 


Uma sucessão de subconjuntos C”. de um espaço 
métrico C diz-se completa, relativamente a um outro sub- 
conjunto S de C, se for possível determinar um inteiro 
positivo N tal que, para n —N, exista pelo menos um 
elemento de C”, na vizinhança : de cada elemento 
de S. 


Os teoremas de convergência podem basear-se 
neste conceito geral de completidade. A discussão da 
convergência torna-se mais simples porém, em regra, se 
o elemento de C”, que aproxima um dado elemento e 
ce S for a sua própria T, “imagem em H.. 


O conceito de completidade que resulta de tomar 
para elementos aproximantes as imagens dos elementos 
ce S é pois uma particularização do conceito geral de 
completidade, à qual pode dar-se o nome de T — com- 
pletidade. 


Uma sucessão de subconjuntos C', de C será 
assim chamada T — completa, relativamente a um outro 
subconjunto S de C, se, dado :, for possível deter- 
minar N tal que, para n > N,a T, — imagem de qual- 
quer elemento e de S esteja contida na vizinhança 
de e, isto é, se; 

d(e,T, (e)) -« : para qualquer e E S (4.1) 

É importante notar que uma sucessão T — com- 
pleta será simplesmente completa, embora o contrário 
não se verifique necessaáriamente., 


Os teoremas de convergência podem ser enun- 
ciados com vários graus de generalidade, correspon- 
dentes a situações mais ou menos complexas. 


Para cada grau aparecem dois teoremas, cada um 
dos quais se relaciona com os dois teoremas de apro- 
ximação da Secção anterior. 


Na presente comunicação, limitar-nos-emos a con- 
siderar um único espaço H, uma única família F e um 
único elemento ec H. 


O primeiro teorema de aproximação dá então ori- 
gem ao seguinte teorema de convergência: 


«Se 
a) os operadores T, forem conformantes e C', CC 
(condição de conformidade), 


b) a sucessão [ C', ' for completa relativamente 
a um subconjunto S de C que contenha s (condição 
de completidade), 


então, (San) converge para 5», 


As condições a) e b) correspondem com efeito às 
condições a) e b) do primeiro teorema de aproximação, 
e não há portanto necessidade de comentários. 

O segundo teorema de convergência, que resulta 
do segundo teorema de aproximação, é o seguinte: 
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«Se 


a) a sucessão 16d for completa relativamente 
a um subconjunto S de C que contenha s, 


b) (S'ank convergir para um elemento de S, 
então [84n) converge para s». 


Também aqui não há necessidade de justificar. 
De acordo com o estabelecido na Secção anterior, 
a condição b) será verificada se a sucessão de ele- 
mentos s,, tais que T, (s,,) =S',n convergir para 
um elemento de S. Ou então, se a sucessão das pro- 
jecções p, dos elementos s',, sobre C convergir para 


um elemento de S. 


n 


5 — Discretização 


O problema da determinação de elementos apro- 
ximantes, simplifica-se extraordinâriamente se o opera- 
dor T for tal que o subespaço H' tenha um número 
finito de dimensões. O problema variacional diz-se então 
discretizado. 

Nesse caso, e se o operador for conformante, tem- 
-se o método de Ritz, cuja convergência pode pois ser 
considerada à luz do primeiro teorema da Secção 
anterior. 

Se o operador não for conformante, tem-se um 
método mais geral que o de Ritz. Trata-se de uma 
situação que aparece frequentemente quando a discre- 
tização se faz pela técnica dos elementos finitos [4, 5] 
que utiliza funções coordenadas definidas seccional- 
mente no domínio a que se refere o problema varia- 
cional. Serve então o segundo teorema, 


6 — Conclusões 


À teoria que acaba de ser apresentada corresponde 
a uma generalização das análises de convergência que, 
para o caso linear, têm sido apresentadas pelo autor em 
vários trabalhos dedicados ao método dos elementos 
finitos [5,6] e à Teoria das Estruturas [7,8] em geral. 

O primeiro teorema de convergência pode consi- 
derar-se clássico, no sentido de que está, de certo 
modo, substancialmente incluído na teoria da conver- 
gência do método de Ritzl3], 

Não assim o segundo teorema, o qual já permitiu 
mostrar que é possível conseguir convergência no mé” 
todo dos elementos finitos mesmo que se viole a con- 
dição da conformidade entre elementos, e explicar a 
razão por que, em certos casos, a convergência não 
chega a obter-se 

A importância deste teorema é porém muito maior 
se se atender ao papel essencial que o conceito de con- 
vergência desempenha na justificação de teorias apro- 
ximadas e dos métodos de geração dessas teorias [7. 8]. 

Assim, na teoria das Estruturas, por exemplo, os 
teoremas de convergência levam a uma justificação dos 
métodos energéticos de estabelecimento das equações 
das teorias das peças laminares e lineares, métodos 
esses que, apesar da sua comodidade, têm sido olhados 
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com desconfiança, em virtude de aparecerem com 
carácter heurístico mais ou menos disfarçado, Mostrar- 
-Se-á num próximo trabalho como o segundo teorema 
de convergência pode ser utilizado na justificação da 
utilização desses métodos energéticos. 


REFERÊNCIAS 


[1] — Vainberg, M. M. — Variational Methods for the 
study of Non — Linear Operators, Holden Day, 
1964. 

[2] — Taylor, A. E. — Introduction to Functional Ana- 
lysis. John Wiley, 1958. 

[3] — Mikhlin, S. G. — Variational Methods of Mathe- 
matical Physics. Pergamon Press, 1963. 

[4] — Oden, J. T. — Finite Element of Non-linear Con 


Arantes e Oliveira, E. R. — «Theoretical Foun- 
dations of the Finite Element Method». Int. J. 
Solids Structures, Vol. 4, No. 10, October 1968, 
pp. 929-952. 

[6] — Arantes e Oliveira, E. R. — «Completeness and 
Convergence in the Finite Element Method». 
Proc. 2nd. Conf. on Matrix Methods in Structu- 
ral Mechanics (AFFDL-TR-68-150),  Wright- 
-Patterson AFB, Ohio, 1968. 

[7] — Arantes e Oliveira, E. R. — «The Convergence 
Therems and their Role in the Theory of Struc- 
tures». Proceedings of the IUTAM Symposium 
on High Speed Computing of Elastic Structures, 
Liêge, 1971. 

[8] — Arantes e Oliveira, E. R. — Lições do Curso 

da NATO sobre «Elementos Finitos na Mecã- 

nica do Contínuo». Edição provisória do LNEC, 


[5] 


tinua. Mc Graw-Hill, 1972. 


Lisboa, 1971. 


INII — Instituto Nacional de Investigação Industrial 
Formação 72/73 


Programas e Datas 


DPE.33 — Psicologia So- 
cial e Técnicas de comu- 
nicação 


12 a 16 de Março 


35 horas 


S1.31 — Matemática e Es- 
tatistica de Gestão 


Início em 12 de Março 
(3 sessões por semana em 
fins de tarde) 


36 horas 


S1.32 — Princípios Gerais 
de Direcção 


12 a 16 de Março 


30 horas 


SEC.31 — A Organização 
de Empresas e de Estalei- 
ros de Construção Civil 


8 a 10 e 12 a 16 de 
Março 


32 horas 


Tema 


A empresa e os problemas 
da comunicação. Contri- 
buição da psicologia so- 
cial. 


As bases matemáticas da 
gestão da empresa, 


Conceitos e noções funda- 
mentais do management. 
Estratégia directiva e or- 
ganização do pessoal diri- 
gente. 


A organização racional da 
indústria da construção ci- 
vil. 


Objectivos 


Mostrar a importância do 
fenómeno relacional nas 
empresas e chamar a aten- 
ção para os problemas 
globais que daí resultam. 


Apresentar os principais 
instrumentos matemáticos 
e estatísticos que intervêm 
mais frequentemente na 
gestão. 


Situar o participante na 
sua organização e na so- 
ciedade. Levá-lo a um me- 
lhor conhecimento das 
funções mais importantes 
da organização e da sua 
integração numa direcção 
eficiente, Desenvolver as 
aptidões de direcção e de 
trabalho em grupo. 


Equacionar os problemas 


que as empresas e os es- 
taleiros de construção ci- 
vil enfrentam no domínio 
da organização científica 
e apontar vias para a sua 
solução. 


Dirigido a: 


Dirigentes e quadros su- 
periores das empresas e 
dos serviços, nomeada- 
mente os responsáveis pe- 
los seus aspectos huma- 
nos. 


Todos os quadros que 
exerçam funções de res- 
ponsabilidade  participan- 
do na preparação de deci- 
sões na sua organização. 


Quadros intermédios na 
hierarquia das empresas 
ou sreviços, quadros com 
a dupla função de contro- 
lar os seus departamentos 
e de coordenar as suas 
actividades com as dos 
outros. 


Empresários, dirigentes e 
quadros técnicos das em- 
presas da construção civil. 


Para informações mais detalhadas sobre estas acções e respectivos programas, os interessados poderão dirigir-se ao 
INII — Divisão de Formação e Aperfeiçoamento, Rua Filipe Folque, 2, 6.º-D. — Lisboa-1 — Telefs. 53 36 20/3/4/5/6. 
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SOBRE A REOLOGIA DE AREIAS 


RESUMO 


A partir da vibração em ressonância de tubos me- 
tálicos cheios com material granular, é possivel deter- 
minar os parâmetros reológicos do enchimento. Um 
ensaio com uma areia de mármore mostrou que se pode 
atribuir a esta areia comportamento visco-elástico (corpo 
de Maxwell). O método parece aplicável a problemas 
de Mecânica dos Solos. 


1 — INTRODUÇÃO 


A investigação do comportamento reológico de 
substâncias granulares foi iniciada com o fim de se 
conhecerem as propriedades do manto terrestre (que se 
supõe formado por rochas peridotíticas) quando sujeito 
a fusão parcial. O método usado interessa porém igual- 
mente para o estudo reológico de qualquer areia. 

As substâncias granulares têm propriedades que 
são, de certo modo, intermédias entre as dum sólido e 
as dum líquido, parecendo razoável admitir que o mo- 
delo reológico correspondente seja (pelo menos em 
primeira aproximação) o corpo visco-elástico de Maxwell. 
O facto terá naturalmente de ser confirmado pela expe- 
riência. 

Uma substância visco-elástica é, como se sabe, 
definida pela relação (veja-se, por exemplo, Jaeger, 
1956, p. 102) 


2nd/dt 


E . (1) 
1+ (rp) d/dt 


onde ,» e , são os tensores distorcionais («deviators») 

respectivamente das tensões e das deformações, t é o 

tempo, 7 a viscosidade e u,, a rigidez de Maxwell. 
Num sólido elástico perfeito seria 


q=2 ue (2) 
sendo » o módulo de rigidez habitual (módulo de 
elasticidade transversal). 


Pode dizer-se portanto que num corpo de Maxwell 
o módulo de rigidez é substituído pelo operador 
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F. MACHADO 
(Dr. em Engenharia Civil, I.S.T.) 


e M. E. CORTE-REAL 
(Aluna de Eng.º Electrotécnica, 1.5.7.) 


ABSTRACT 


A method (based on resonant vibration) was deve” 
loped for measuring the rheological constants of sandy 
material contained in small metallic pipes. A test on a 
marble sand showed that a viscoelastic behaviour 
(Maxwell body) can be ascribed to this material. The 
method is expected to be useful for some problems of 
Soil Mechanics. 


«d/dt 
 1+(n/uy)d/dt 


(3) 


Por sua vez, para um líquido viscoso newtoniano 
seria 


c=2nde/dt (4) 
o que é equivalente agora ao operador 
un d/dt (5) 


Numa vibração em que ; e . são ambos pro- 


mit! 


porcionais à parte real de e x “(sendo + O período), 


virá 


d/dt=2 = i/1 (6) 


e as expressões (3) e (5) passam a escrever-se, res- 
pectivamente, 


2min/z 


1+25in/um* o 


LM Ha T 
= — Sema (15 | (7) 
| + vu? 2/472%2 2%7n 
p= 2min/r (8) 


As equações (2), (7) e (8) podem considerar-se 
casos particulares da expressão geral 


n=puo(1+i/0) (9) 
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onde », é o módulo de rigidez (correspondente ao 
período considerado) e 1/Q é um coeficiente de amor- 
tecimento (OQ é o chamado factor de qualidade em 
ressonância). 

Para um sólido perfeito será 


Lo 1 1/Q = 0] (10) 


para uma substância visco-elástica de Maxwell 


u | Una 7 
qa E. ds (11) 
o | um? 72/4752 %2 Q 277 
e para um líquido newtoniano 
“4 =0 1/Q= o mas hs I0=2mn/)/% (12) 


É interessante notar que, na substância visco-elás: 
tica, quando o período de vibração for muito pequeno 
1/Q é vizinho de zero e a substância comporta-se apro- 
ximadamente como sólido elástico de rigidez vu. 
Quando, pelo contrário, o período for muito grande u, 
aproxima-se de zero e 1/Q de infinito, mas o produto 
to /Q conserva o valor finito 2 = 7 | 7, como sucede 
nos líquidos newtonianos. As substâncias visco-elásticas 
podem assim considerar-se corpos de tipo mais geral 
que incluem como casos extremos os sólidos elásticos 
(perfeitos) e os líquidos viscosos (newtonianos). 


| LIQUIDO 


(QUASE 
NEWTONIANO ) 


H/Q 


SOLIDO 


(QUASE ELÁSTICO) 


0 : /T 


Fig. 1 —- Variação teórica de »/Q em função de 1/7, para corpo 
de Maxwell 


As considerações anteriores mostram haver inte- 
resse em representar a relação p. / OQ como função de 
1/= (Fig. 1). O máximo da curva estabelece, de certo 
modo, a separação entre o comportamento de líquido 
viscoso e o de sólido elástico. O período 7. — 27=7%7/tm 
correspondente ao máximo coincide com o chamado 
tempo de relaxação. 

Atendendo a (11), a grandeza O/p. + é linear em 
1/72, o que facilita o ajustamento pelos mínimos qua- 
drados, conforme se verá adiante (cf. Fig. 5). Na repre- 
sentação gráfica correspondente a ordenada na origem 
é 1/2 77 e o coeficiente angular da recta é 2 xr7/py?. 

Para esforços moderados, as imperfeições de elas- 
ticidade que afectam o módulo |. , não têm influência 
no módulo de compressão y (módulo de elasticidade 
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volumétrica dos sólidos elásticos). Nos problemas de 
vibrações, as soluções obtidas para um sólido perfeito 
(quando expressas nos módulos y e ») podem esten- 
der-se aos semi-sólidos fazendo a simples substituição 
indicada em (9). O método vai ser utilizado na secção 
seguinte. 


2. TEORIA DA VIBRAÇÃO DE TUBOS 


O dispositivo experimental para estudar a resso- 
nância de barras vai indicado na Fig, 2 (veja-se, por 


TRANSDUTORES 
a" 


AMPLIFICADOR 
DE VOLTAGEM 


OSCIL ADOR 
ELECTRICO 


AMOSTRA 


OSCILOSCÓPIO 
DIGITAL CATÓDICO 


FREQUENCIMETRO 


Fig. 2 —- Esquema do equipamento usado nas medidas reológicas 


exemplo, Spinner & Tefft, 1961). Uma oscilação eléc- 
trica depois de amplificada é convertida em vibração e 
aplicada à amostra através de um dos fios de suspensão. 
O outro fio está ligado a um detector que transforma 
agora a vibração em oscilação eléctrica, a qual depois 
de amplificada vai aparecer no «écran» dum oscilos- 
cópio catódico. 

Enrolando os dois fios em sentidos contrários, ex- 
citam-se não só vibrações forçadas em flexão, mas tam- 
bém em torção. Fazendo variar a frequência do oscilador, 
obtêm-se as várias ressonâncias (que se reconhecem por 
máximos da amplitude na imagem do osciloscópio). Es- 
tas frequências de ressonância são, como se sabe, prã- 
ticamente iguais às frequências de vibração livre da 
barra. 

A fim de aumentar o rigor, a frequência ou o pe- 
riodo da ocilação excitadora é medido num frequenci- 
metro digital (Fig. 2). 

No presente estudo há interesse em utilizar tubos 
de secção circular, À vibração transversal livre dum tubo, 
infinitamente estreito, é descrita pela equação diferencial 
(cf. Timoshenko, 1955, p. 324 e seg.) 


4 2 
) 

je = — po (13) 
dx* dt? 


onde y é o desvio vertical à distância x de um dos ex- 
tremos, t é o tempo, : a densidade, À a área da secção, 
E o módulo de Young e | o momento de inércia da sec- 
ção em relação a um diâmetro. Sendo, respectivamente, 
de ed, os diâmetros exterior e interior dos tubos, virá 


A=z(dç2-d,2)/4e1-z(d4-d,4)/64 
(14) 
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Se nos extremos do tubo não actuarem forcas exte- 
riores, as condições nos limites são 


d?yldx2=o para x=0 e x=1 (15) 
d3y'dx3 ='o para x=0 e x=1 (16) 
onde | é o comprimento do tubo. 
O período da vibração é 
21 |? / - 
- is Pre de (17) 
32 IE 
sendo 5 uma solução de 
cos 3 cosh 5 — 1 (18) 


Para o modo fundamental é 5 = 4,730 e para a 
2.º harmónica 5 = 7,853. 

À equação (17) permite calcular E. Para um tubo 
real (de secção não nula) será 


a SE (19) 


sendo C um factor de correcção que se pode representar, 
com suficiente aproximação, pela fórmula empírica 


C=1+k/12 (20) 


onde k é uma constante (que depende do material do 
tubo e da secção transversal). 

Para tubos do mesmo material e com a mesma 
secção C/E é uma função linear de 1/12, o que permite 
determinar cômodamente o módulo E. 

Se o tubo estiver cheio com uma substância de 
viscosidade suficientemente elevada para vibrar como 
um sólido, será 


onde agora o indice t se refere ao tubo e o índice a 
à amostra que constitui o enchimento. 

Para o tubo cheio pode determinar-se um módulo 
equivalente, usando a equação (19) onde se fará 


l 


A=A +Api=h + 00= A à +AstO) 
(andas + ai 


sendo evidentemente 


IEL E, +I,E, ou dE -(dç4— dE + dE, 


(23) 
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O módulo E, pode obter-se desta equação. O valor 
só é, porém, verdadeiro se forem iguais os coeficientes 
de Poisson dos tubos e do enchimento. No caso contrá- 
rio, o valor determinado é aparente (sendo representado, 
no que se segue, por E, * ). 

Se forem o es, as tensões radial e longitudinal 
no material do enchimento e :, e :, as dilatações 
lineares correspondentes, o equilíbrio elástico implica 
o seguinte sistema de equações (Timoshenko, 1940, 


p. 62; 1941, p. 241) 


l 
Ee E PD == 90) 0, == 946 


a 


Z. 


—. 


' a E, do? — d;? 
(24) 
onde os vv são coeficientes de Poisson, 
Eliminando G er obtemos 
EAtl—A+2myy! 
g cá rs SR A Rd £ =— Es E, 
* (Ao u+ro * 
donde se deduz (por ser 2 p, E /(1:7,)) 
oo GÊ O Bo ta (26) 
» 2 [8º =] + 27)pa| 
com 
=" (1+Bv,/n)/ (1 + B) (27) 
8 


(28) 


Na prática, = “ é, em regra, uma boa apro- 


4 t 
ximação. 


Para as vibrações em torção teremos anâálogamente 
(Timoshenko, 1955, p. 318 e seg.) 


ls - (29) 


sendo agora «w o ângulo de rotação à distância x de um 
dos extremos, » o módulo de rigidez (ou de elastici- 
dade transversal) e | o momento de inércia polar (em 
relação ao centro da secção) que é 


U=21== (d(4— d;4)/32 (30) 
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Para extremos livres, as condições nos limites são 

doldx=o para x=o e x =] (31) 
O período desta vibração é 

ad elu donde u — 4:12,52 (32) 


Para o tubo cheio vem agora 


9º? CS (a) 
Von + I,p (1 Sr lct)—rs 
gt ata) dx? Lo ata” sp 
donde se deduz, com | 1,“ I,, 


Lp -Lu+Hl,t, OU du=(dç*—d*)pr+ dp, 


(34) 


sendo wu a rigidez equivalente do tubo cheio que se 
obtém da equação (32) fazendo 


(35) 


A equação (34) permite naturalmente determinar 

u,- Conhecidos E, * e w, obtém-se, pela equação (26), 

o coeficiente de Poisson »,, e depois, utilizando as 

equações conhecidas, o módulo de Young (verdadeiro) 
E, e o módulo de compressão , 

E. 24 


a 


(144) (36) 


a 


4 =E/(3-64,) (37) 


Para materiais granulares há interesse em usar 
pequenas tampas nos extremos dos tubos, 

Por analogia com barras com discos nos extremos 
(cf. Timoshenko, 1955, p. 319 e seg.) os períodos =, 
ou 7,, dos tubos sem tampas, respectivamente vazios ou 
cheios, estão relacionados com os correspondentes perío- 
dos dos tubos com tampas ( 7, ou ) pelas expres- 
sões aproximadas 


“ta 


<.2 SI TA/A, 8) 


t 


(38) 
Tg — Tuaé (1 Tel A, 6 + As 84) 
para as vibrações em flexão; e 
2 Tu2 / (1 + «/ (1, 5) | 

(39) 
Ta “tra? / [1+: /L, RE sl ta 11 


para as vibrações em torção. 

As constantes : ou : (para cada tubo) são gran- 
dezas pequenas desde que as tampas sejam ligeiras. 
É evidentemente possível eliminar : e : nas equações 
(38) e (39). 
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O resultado pode exprimirse dizendo que os valo- 
res equivalentes E/C, para os tubos cheios (sem tam- 
pas), podem ser obtidos da equação (19) usando o pe- 
riodo correspondente aos tubos cheios e com tampas, 
desde que se substitua a densidade pelo valor 


d,º Tu? d,? 
a ade De RT (so) 


em vez do valor dado por (22). 

Anãlogamente, a rigidez » do tubo cheio (sem 
tampas) pode ser calculada por (32) com o período 
medido com tampas, fazendo, em vez de (35), a subs- 
tituição 


d 4 - d.4 
ç (1 7, — ç + —— e. (41) 


Para estudar o amortecimento, a vibração dos tubos 
pode representar-se por uma elongação proporcional à 
parte real de 


e 2 RINIS (42) 
O efeito de imperfeições de elasticidade obtém-se 


fácilmente substituindo em (17) ou (32), e depois em 
(42), 


E = E 11 + 1/0.) (43) 
ou 
nus 1 +i/0) (44) 
sendo 
] 2 (1 + vu) | 
= 45 
Q, 3 Q io 


que se deduz de E=9,u/(317 u), supondo 1/0 
suficientemente pequeno para, depois da substituição 
(44), se poderem desprezar os respectivos quadrados 
(7 , como se disse, não é afectado pelas imperfeições 
de elasticidade). 

O resultado é o aparecimento de factores de amor- 


tocimento 6 7!/700t ou e TO (To corresponde a 


E, ou w., respectivamente). Numa oscilação forçada 
com estes amortecimentos a amplitude da vibração é 
proporcional a (cf. Timoshenko, 1955, p. 78) 


| | | | ) ] 1. T=H 

e fe (46) 
! -2 To? -2 To? Q, 2 
ou a uma expressão análoga em 1/0, 7, sendo agora 
o período de ressonância e = o periodo da excitação. 


A amplitude é reduzida a 1/y 2 do máximo para 
valores de 7 (de ambos os lados de To ) dados apro- 
ximadamente por 


Qlr=— milto = 1/0 (47) 
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onde se admitiu que + é vizinho de , ; há evidente- 
mente uma expressão análoga para 1/Q (vibração em 
torção). 

isto dá um método para determinar 1/0, ou 1/0, 
reduzindo de 1/1, 2 a amplitude exibida no osciloscópio 
(para ambos os lados de cada ressonância) e medindo 
os períodos correspondentes. 

Substituindo agora em (23) e (34) expressões aná- 
logas a (43) e (44) para os tubos e para o enchimento, 
obtém-se (igualando as partes imaginárias) 


dçs E, 1Q (ds 68) Es: /Q Er E d;* Esa Ora (48) 


E 


des Lo /Q (dç4 e d; +) ot Q d;* Poa Q, (49) 
equações que permitem determinar 1/Q0,, e 1/Q, co- 
nhecidas as restantes grandezas. 

O efeito da suspensão que tanto afecta 1/0, com 
1/Q,, (e é proporcional, em cada caso, a T/IE, ou a 
1/1, E,,) fica eliminado na subtracção. O mesmo acon- 
tece em relação a 1/0 e 1/0,. 

Na prática, para as vibrações em flexão, determi- 
nam-se (em vez de E, /Q,,) valores aparentes E, * / 


/Q,,* que se podem transformar pela expressão apro 
ximada 
Boa! 0, [13/89/01 + 1,)] Et / Ogg (50) 


Isto resulta de ser 

Eq tas) ra + (4/3) 1 + 7)pa (51) 
onde se pode fazer y — » e desprezar a parte imagi- 
nária de » (o que é uma razoável aproximação porque 


o amortecimento dos tubos é, em geral, pequeno, com- 
parado com o do enchimento). 


3. APLICAÇÃO AO ESTUDO DUMA AREIA 


Utilizou-se areia artificial preparada com mármore 
de Estremoz; aproveitou-se a fracção que passava no 
peneiro n.º 10 e era retida no n.º 35. A densidade deste 
mármore é 2,695 g/cmS3. 

Usaram-se quatro tubos de aço inoxidáxel com as 
dimensões e densidades indicadas na Tabela 1. As tam- 
pas eram de plástico e encaixavam bastante justas. 

Para cada tubo (sem tampas) mediram-se os pe- 
ríodos das ressonâncias fundamental e 2.º harmónica 
em flexão, e fundamental em torção, Com esses perio- 
dos calcularam-se, pelas equações (19) e (32), os valo- 
res E/C e u» que constam da Tabela 2. 

Na Fig. 3 representaram-se as grandezas C/E em 
função de 1/12; como era de esperar, os pontos definem 
duas rectas cujo ajustamento pelos mínimos quadrados 
permitiu determinar o módulo E. Para » adoptou-se a 
média aritmética dos valores da Tabela 2. Os resultados 
são 

E = 1,956 + 0,003 Mb 


u = 0,7526 + 0,0003 Mb 
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a que corresponde (usando uma equação análoga a (36)) 
o coeficiente de Poisson 


v = 0,299 + 0,003 


05” 


TUBOS VAZIOS | | 


07+ a 4 4 — 
| 


2* HARMÔNICA 


= | 
E | 
> 06 a 
uu) 
o 
Q FUNDAMENTAL | 
os 
| 
do 0005 001 0015 002 0025 
1/0 (em?) 
Fig. 3 —- Diagrama de C/E (em função de 1/12) para os tubos 
vazios 


Os amortecimentos foram medidos para os tubos 
com tampas. Os valores apresentavam alguma dispersão 
mas foram regularizados (usando a teoria desenvolvida 
noutra parte; Machado, artigo no prelo); isto tem, de 
resto pouca influência nos resultados finais. Os valores 
adoptados indicam-se na Tabela 3. 

Determinaram-se também os períodos de resso- 
nância dos tubos vazios com tampas, para poderem usar- 
-se as equações (40) e (41). A relação 7, /7, (para as 
tampas utilizadas) variava entre 1,041 e 1,066. 

Os tubos encheram-se depois com a areia, a que 
se deu uma ligeira compactação. A densidade aparente 
obtida foi, em média, de 1,33 g/cm3, o que corresponde 
a uma porosidade de cerca de 51 %. 

Determinaram-se então as frequências de resso- 
nância dos tubos cheios, as quais pelas equações (19), 
(32), (40) e (41) conduzem aos valores E/C e u indi- 
cados na Tabela 4, 

Na Fig. 4 apresentam-se os valores C/E em função 


09: r ,-—- , T 


TUBOS CHEIOS 


08 


* HARMÔNICA 


C/E (Mb!) 


0.6: 
“FUNDAMENTAL 


054 0005 001 0015 002 0025 


VC (em?) 


Fig. 4 —- Diagrama de C/E (em função de 1/12) para os tubos 
cheios 
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de 1/12. Procedendo como para os tubos vazios, ob- 
têm-se os módulos equivalentes 


E — 1,621 + 0,006 Mb 
u = 0,579 + 0,002 Mb 
e pelas equações (23) e (34) deduzem-se para a areia 
E* — 0,746 Mb 
4 — 0,123 Mb 


Este é um valor aparente do módulo de Young. 
Pela equação (26) pode determinar-se o coeficiente de 
Poisson 


1. = 0,358 


ua 
e pelas equações (36) e (37) 


E, — 0,335 Mb 


74 = 0,394 Mb 


Conforme veremos, a areia tem comportamento que 
se pode considerar visco-elástico e por isso », varia com 
o período de vibração (de acordo com a equação (11), 
o mesmo sucedendo a E, e »,. Os valores apresen- 
tados representam, de certo modo, médias no intervalo 
dos períodos de ressonância utilizados. O módulo 7, 
calculado pode, porém, ser vizinho do valor verdadeiro 
(se E,* e », corresponderem sensivelmente a um mes: 
mo período de vibração), 


Os amortecimentos dos tubos cheios (e os respec- 
tivos períodos) vão indicados na Tabela 5. Não foi pos- 
sível medir 1 /Q (ou 1/0, ) em muitas das ressonân- 
cias, ou por ser muito pequena a amplitude da vibração, 
ou por haver sobreposição doutras ressonâncias anóma- 
las (que resultam certamente da heterogeneidade da 
areia e eventualmente de efeitos da suspensão). 

Dos valores medidos obtiveram-se pelas equações 
(48), (49) e (50) os produtos toa/Q, que permitiram 
representar na Fig. 5 Q/v 7 em função de 1/72, Os 


03 


ot 


Q/ut (10ºMb's”) 


1/7? (10º87*) 


Fig. 5 — Diagrama de QO/ 47 (em função de 1/ 7 2 | para a areia 


122 


pontos definem razoávelmente uma recta de inclinação 
finita e não nula, o que mostra que o comportamento 
reológico da areia pode ser descrito, dentro da aproxi- 
mação adoptada, pelo modelo visco-elástico de Maxwell. 

Um ajustamento pelos mínimos quadrados permite 
calcular a ordenada na origem e a inclinação, e a partir 
destes valores os parâmetros de Maxwell (cf. Fig. 5) 


um= 0,16 + 0,03 Mb 


x 
| 


(2,7 + 0,7) X 108 poises 


O tempo de relaxação correspondente é da ordem 
de 10 —4s apenas; por isso, mesmo para solicitações 
de duração relativamente curta, as areias escoam como 
líquidos viscosos. 

Com os parâmetros obtidos é possível determinar, 
para os vários períodos de vibração, valores de p,, dados 
por (11). Na Fig. 6 a curva teórica de variação de p, 


020 + |I | - - ' 1 
| | 


nd O] 
e 
4 


H (Mb) 


CURVA TEORICA 
01 


25 l so 75 100 


T(us) 


Fig. 6 -- Diagrama de |. em função de 7, para a areia. À curva 
representa a variação teórica (a curva tracejada corresponderia a 
um = 0,165 Mb, “ =3,0 X 10º poises) 


é comparada com os valores » | obtidos directamente da 
Tabela 4, aplicando a cada medida de wu a equação 
(34). Os pontos acompanham razoávelmente a curva 
(que é de resto, afectada pela incerteza de ny, e 1),0 
que está também de acordo com o modelo visco-elástico 
considerado. 


4. CONCLUSÕES 


Conforme era certamente de esperar, o comporta- 
mento da areia ensaiada pode ser descrito, com razoável 
aproximação, pelo modelo visco-elástico de Maxwell. Os 
respectivos parâmetros reológicos são 


0,39 Mb 


| 


0,16 Mb 


Mm 


«=> 2,7 X 108 poises 
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O método usado parece suficientemente simples 
para ter interesse prático em problemas de Mecânica dos 
Solos. 
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TABELA 1 


Dimensões e densidades dos tubos 


Comprimento Diâm. ext. Diâm. int. Densidade 
(cm) (em) (cm) (g/cm3) 
Tubo A 7,0698 0,8031 0,5821 7,946 
Tubo B 8,1582 0,8033 0,5824 7,970 
Tubo C 10,0024 0,8025 0,5823 7,952 
Tubo D 14,1190 0,8020 0,5811 7,923 
TABELA 2 
Constantes elásticas dos tubos (vazios) 
E/C, fund. E/C, 2.º harm. “, fund. 
(Mb) (Mb) (Mb) 
Tubo A 1,748 1,431 0,7522 
Tubo B 1,793 1,535 0,7532 
Tubo C 1,851 1,654 0,7519 
Tubo D 1,899 1,789 0,7531 
TABELA 3 
Amortecimentos 1/0 dos tubos vazios com tampas 
1/Q ;, fund, 1/0, 2º harm. 1/0, fund. 
(10 - 3) (10 - 3) (10 - 3) 
Tubo A 1,109 1,004 1,166 
Tubo B 1,171 1,008 1,169 
Tubo C 1,320 1,021 1,174 
Tubo D 1,910 1,070 1,186 
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Constantes elásticas (equivalentes) dos tubos cheios de areia 


TABELA 4 


E/C, fund. E/C, 2.º harm. u, fund, 

(Mb) (Mb) (Mb) 
Tubo A 1,537 1,227 0,586 
Tubo B 1,565 1,308 0,579 
Tubo C 1,584 1,386 0,576 
Tubo D 1,597 1,507 0,576 

TABELA 5 

Amortecimentos 1/Q dos tubos cheios de areia 

Período 1/0, 1/0 

(us) (10 - 2) (10 - 2) 
Tubo B 51,8 — 3,21 
Tubo C 67,1 — 3,41 
Tubo C 96,0 3,04 — 
Tubo A 129,7 1,92 — 
Tubo C 250,1 2,25 — 
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QUANTO MAIS 
FORMOS CONHECIDOS 


Milo oERENDO 
APREGIADDO 


N 
ana na a A A AT À 


5 »” 
AU cs as 
“MM YWVv7 Es 


4 HD TT —e—eeoo 


“ 
Ed 
sa” 


| E 


o no aa e 
= "UMA .- 
AEB & cr» BAIA c A A + 


e a e e O O e it e CS SSD e e e —————e +. 


O certo é que contribuimos para que um sem 
número de produtos e serviços básicos 
vos cheguem às mãos eficientemente. 


A luz eléctrica, o gás, a gasolina, os óleos, 
o aço, O cimento, o papel e até o próprio avião 
estão de certa forma ligados à nossa 
empresa. E que ajudamos a construir e a criar 
indústrias e serviços de base. 


Para o grande público talvez sejamos anónimos, 
mas para o desenvolvimento industrial 
do país somos indispensáveis. 
E porque servindo eficientemente o país 
também o servimos a si, 
gostariamos que nos conhecesse melhor. 


Sabia por exemplo que construímos 
e montamos as blindagens do túnel e da 
e aT=aglia lho /=Kcro |U |||og/oN-=E-N oo pfoBir-hjo go -[o/- 

de Vilarinho das Furnas, 
oops ooo lh rojr-| No |-0M=]0/0 Mb voa]=|[-[o[=[=Tá 
E doravante sempre que passar por uma 
barragem ou uma central eléctrica... 
talvez se lembre do nosso nome. 


Sepsa 


SOCIEDADE DE CONSTRUÇÕES ELECTROMECANICAS, S.A.R.L. 


TELEFS. 951616 (10 LINHAS) 
APARTADO S-LECA DO BAILIO - TELEX: SEPSA P 2616 


| 
servindo eficientemente o país, também o servimos a si... 
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MICROSCOPIA QUANTITATIVA DE CLÍNQUERES E CIMENTOS PORTLAND 


RESUMO 


A microscopia óptica com luz reflectida apresenta 
grande interesse no estudo do clinquer e do cimento 
portland pelo seu duplo aspecto: 

— Qualitativo, facultando a observação directa da 

estrutura e textura dos grãos de clinquer. 

— Quantitativo, possibilitando a determinação di- 

recta das proporções dos seus constituintes. 

No presente trabalho, após breve revisão dos aspec- 
tos fundamentais respeitantes à composição do clinquer 
e análise microscópica qualitativa, descreve-se a téc- 
nica de determinação quantitativa das fases pelo método 
de contagem de pontos com uma ocular integradora 
de reticulo de pontos. 


1 — INTRODUÇÃO 


A microscopia óptica pode considerar-se já como 
uma técnica tradicional na análise de cimentos. Le Cha- 
telier (1882), pela primeira vez, observou o silicato tri- 
cálcico numa lâmina delgada de clinquer com um 
microscópio polarizante e depois Tornebohm (1897) 
identificou e descreveu os constituintes mineralógicos 
principais do clínquer. O método por ataque selectivo 
em secções polidas de grão de clínquer só foi aper- 
feiçoado mais tarde por Tavasci (1934) e Insley (1936). 
Desde então, o microscópio metalográfico tem sido 
utilizado com o máximo interesse e êxito na investiga- 
ção e no controle industrial do cimento portland, 

O exame microscópico do clínquer informa sobre 
a espécie, formação e repartição das fases e fornece 
indicações importantes sobre os fenómenos que podem 
ocorrer no decurso da fabricação do clínquer e prepa- 
ração do cimento. O emprego duma platina ou ocular 
integradoras faculta ainda, por um processo de conta- 
gem mineralógica, uma determinação quantitativa directa 
dos constituintes integrantes do material. 

Foi já objecto dum relatório do LNEC um estudo 
microscópico qualitativo de clinqueres e cimentos nacio- 
nais [3]. Por isso, no presente trabalho, apenas se 
passarão em revista alguns aspectos fundamentais da 
microscopia qualitativa, descrevendo-se mais em por- 
menor o processo de determinação quantitativa por con- 
tagem de pontos com ocular integradora, recentemente 
adquirida. 
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SYNOPSIS 


The reflected light optical microscopy is of great 
interest for the study of clinker and portland cement, 
both from the qualitative viewpoint since it allows the 
examination of the structure and texture of the clinker 
grains, and from the quantitative viewpoint as it is pos- 
sible to determine the proportions of its constituents. 

In the present work, after a brief review of the fun- 
damental aspects related to clinker composition and to 
qualitative microscopic analysis, a description is made 
of the technique of quantitative determination of phases 
by the point counting method using a integrating eyepiece 
with point grid. 


2 — CONSTITUIÇÃO MINERALÓGICA DO CLINQUER 
DE CIMENTO PORTLAND 


O clínquer pode considerar-se como uma rocha 
artificial cuja formação se deve a um conjunto de 
reacções em estado sólido com fusão parcial dos com- 
ponentes. O cimento portland é um ligante hidráulico 
resultante da mistura de clíinquer e gesso finamente 
moídos. A composição das fases do clínquer, a estru- 
tura dependente da velocidade de arrefecimento e a 
finura, entre outras propriedades, definem em grande 
medida o comportamento do cimento. 


Os constituintes mineralógicos, ou fases, normais 
no clínquer podem incluir-se em três grupos distintos: 


1 — Silicatos, em cristais bem formados, que não 
chegaram a fundir durante o período de clin- 
querização. 


2 — Fase intersticial, correspondente à fase fun- 
dida à temperatura de clinquerização, no seio 
da qual cristalizaram os silicatos e que por- 
tanto solidificou posteriormente âqueles. Esta 
fase aglutina os cristais dos silicatos. 


3 — Óxidos, grupo menor que forma inclusões nos 
grupos anteriores 


O 1.º e 2.º grupos incluem os chamados cons- 


tituintes principais do clínquer: alite, belite, aluminato 
e ferrite. 
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2.1 Silicatos 


Neste grupo estão a alite e a belite, as fases exis- 
tentes em maior proporção. 

Por «alite» designam-se as diferentes modifica- 
ções do silicato tricálcico estabilizadas no clinquer im 
dustrial por formação de soluções sólidas com inclu- 
são de pequenas quantidades de óxidos estranhos (MgO, 
Al, O, Fe, 04, TiO,, P, Os, SO, ...). Conhe- 
cem-se actualmente 6 formas alotrópicas do C,S (1 
trigonal, 2 monoclínicas e 3 triclínicas) identificáveis 
por difracção de raios X e análise térmica diferencial, 
de estruturas cristalinas muito semelhantes, vizinhas 
da forma trigonal estável a alta temperatura [4, 5, 18]. 

A alite corresponde ao C,S do processo de cál- 
culo da composição potencial. Ocorre normalmente em 
cristais idiomórficos, pseudo-hexagonais, de bordos por 
vezes corroídos, denunciando reacção com a fase li- 
quida, mais ácida, durante o arrefecimento do clinquer, 


frequentemente com inclusões de belite (fig. 1). A alite 
só é estável acima de cerca de 1200C e por isso, com 
um arrefecimento muito lento, pode cecompor-se em beli- 
te secundária e cal livre secundária (C,S » C,S+ CO), 
A velocidade de decomposição é lenta e portanto 
aquela não se dá, em geral, durante o processo de 
arrefecimento industrial. Não obstante, a inclusão de 
certos iões estranhos, como o Fe2+, pode acelerar 
a reacção e reconhecem-se ao miscroscópio associações 
de belite e cal pseudomórficas com a alite, muitas vezes 
com restos de alite não decomposta (fig. 2). 

Belite é o nome dado ao silicato dicálcico presente 
no clinquer, Conhecem-se 4 variedades alotrópicas do 
C,S, 7, 7, B e y, estáveis em domínios de tem- 
peratura determinados, as quais apresentam estruturas 
cristalinas nitidamente diferentes, Podem representar-se 
esquemáticamente as transformações do C S do se- 
guinte modo [6]: 


680º - 630º 
L eee H 
[cristais 690º + 20 
grandes) 
1425º210 
À — Cu 11609210 
a <500º 
6307-5209 
— ca e 
G, — BL 
(cristais 690º220 


pequenos) 


- 850º- 900º 


Mistura 


estequiométrica 


| 780º - BCOS ba 
e Rd 


Fig. 1 — Ataque nital 


400 x 
Cristais de alite de bordos ligeiramente corroídos e com inclu- 
sões de belite 
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Fig. 2 — Ataque nital 400 x 
Alite decomposta. Reconhecem-se os contornos primitivos e no- 
tam-se ainda restos de alite não decomposta 
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A belite do clinquer corresponde a um C,S-5 de 
rede cristalina ligeiramente modificada pela inclusão 
de óxidos em solução sólida (MgO, BaO, Mn, O,, 
Cr, 0,, K,O, Na,0, SiO,, SO, ...). 

O aspecto dos cristais varia bastante com a com 
posição e as condições de cozedura. Normalmente ocorre 
em grãos arredondados com estriação cruzada caracte- 


A magnésia participa na formação da fase fundida 
e entra em solução sólida nas diferentes fases descritas, 
mas para teores em MgO superiores a cerca de 2% 
cristaliza independentemente como periclase. 

A cal livre (primária) resulta quase sempre de defi- 
ciências de homogeneização do cru e por vezes duma 
sobredosagem em cal. 


rística (fig. 3 e 4). 


400 x 
Belite estriada denotando tendência para a forma «dedo de luva» 


Fig. 3 —- Ataque nital 


2.2 — Fase intersticial 


Engloba as duas fases cristalinas designadas cor- 
rentemente por aluminato tricálcico (C , A) e aluminofer- 
rato tetrálcico, ou simplesmente ferrite (C, AF). Na rea- 
lidada cada um destes constituintes pode diferir da com- 
posição estequiométrica, quer por proporções variá- 
veis de A e F, quer por introdução de impurezas em 
solução sólida [4, 5, 18]. 

O aluminato tricálcico, quando bem diferenciado 
na fase intersticial, apresenta-se em cristais irregulares 
ou pseudocúbicos. Cristais idiomórficos alongados, in- 
correctamente designados por «fase prismática», são 
atribuídos a uma fase alcalina de fórmula NC, As. 

A fase ferrite corresponde a um termo duma solu- 
ção sólida de fórmula geral C, A x F,. x: À com- 
posição desta fase nos clinqueres depende do trata- 
mento térmico e da relação A/F. A forma mais está- 
vel parece ser o C, AF. Esta fase distingue-se bem 
pelo seu elevado poder reflector (fase intersticial clara) 
(fig. 5). 


2.3 — Óxidos 


No clinquer podem aparecer, como constituintes 
acessórios, periclase e cal livre (fig. 5 e 6). 


TÉCNICA 418 


Fig. 4 —- Ataque nital 
Cristais lisos e irregulares de belite 


3 — PREPARAÇÃO DAS AMOSTRAS PARA ANÁLISE 
MICROSCÓPICA 


Uma escolha criteriosa da amostra tem natural- 
mente uma importância primordial dado que esta deve 
ser representativa do material a ensaiar. 


Realiza-se uma primeira amostragem, como é indi- 
cado para a análise química, sobre a amostra de mate- 
rial recebido para ensaio a fim de obter porções mais 
pequenas destinadas a observação microscópica. 


A amostra de cimento assim obtida é utilizada tal 
qual para a confecção dos provetes, A amostra de clín- 
quer é primeiro britada muito cuidadosamente até um 
tamanho de grão inferior a 2,38 mm (fracção que passa 
no peneiro n.º 8 ASTM). Durante a britagem deve 
evitar-se considerável formação de pó, fazendo penei- 
rações frequentes. 


Para a preparação dos provetes, o material deve ser 
transformado, por impregnação num meio próprio, num 
produto compacto e duro que não se desagregue durante 
as operações de desbaste e polimento, indispensáveis 
para obter secções bem polidas aptas ao exame micros- 
cópico em luz reflectida. 


Pode encontrar-se uma descrição pormenorizada 
dos diferentes processos de impregnação, polimento 
e ataque das superfícies na bibliografia citada [7, 8, 9]. 
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400 x 


Fig. 5 —- Ataque nital 400 > 
Zona alítica com cal livre, ferrite e aluminato «prismático» 


Em linhas gerais o modo operatório adoptado é o 
seguinte: 

1 — Impregnação de pequena quantidade da amos- 
tra do material em resina sintética Stratyl A 116 e cata- 
lisador Butanox, em pequenos tubos cilíndricos de 
vidro ( -—- 2 cm. (5). Para o cliquer, esta operação é rea- 
lizada num exsicador de vácuo para auxiliar a liber- 
tação do ar ocluso nos poros dos grãos, Os tubos são 
deixados na estufa a 40 “C durante 24 h. 

2 — O vidro dos tubos é partido e o fundo dos 
provetes é sujeito a desbaste em lixas de esmeril de 
granulometrias sucessivamente decrescentes: 3; 2; 71; 
0/0; 3/0; 4/0, até se obter uma superfície plana e prã- 
ticamente isenta de riscos. 

3 — Polimento preliminar com pasta de diamante 
«Hyprez» de graus 6-, 1º, 1/4-0S sobre panos de polir 
próprios adaptados ao disco duma polidora mecânica, 
empregando como lubrificante o fluído «Hyprez», 
tipo OS. Este polimento tem como finalidade conseguir 
uma superfície polida por corte dos cristais de clínquer 
sem se originarem deformações dos mesmos. 


128 


Fig. 6 —- Ataque água 315 > 


Cal livre 


4 — Polimento final obre um pano de feltro embe- 
bido com suspensão de alumina z em álcool absoluto 
e limpeza final sobre um pano de feltro seco, colocado 
sobre uma placa de vidro. Por vezes, principalmente 
no caso do clíinquer, mais poroso, o provete é lavado 
em álcool absoluto, cerca de 2 minutos, num aparelho 
Ultrasonic mod. 100 L & R, para eliminar o mais 
possível os resíduos do polimento introduzidos nos 
poros. 

5 — Ataque da superfície polida do provete com 
reagentes químicos seleccionados para pôr em evi- 
dência as características estruturais. As diferentes fases 
tornam-se visíveis por ataque selectivo ou pela cor que 
adquirem, dado terem composições químicas diferentes, 
o que motiva velocidades de ataque diversas. A escolha 
do reagente é feita de acordo com o objectivo preten- 
dido. Para análise quantitativa mostrou-se mais vanta- 
joso um ataque com água destilada (4-8 seg.) seguido 
de ataque com nital (3-6 seg.). No quadro 1 estão reu- 
nidos os reagentes mais usuais. 
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QUADRO 1 


Reagentes de ataque mais usuais em análise microscópica de clinqueres 


SS 


Reagentes Concentração Tempos de Acção do ataque 
ataque 
Água = 2-3 seg. Ataca preferencialmente a cal livre 
destilada 
» 3 seg. Começa a atacar as demais fases 
Etanol- Mistura em 1-2 min. Ataca a cal livre 
-água partes iguais 
Nital Solução HNO, 2 - 6 seg. Reagente universal. Ataca todos os constituintes 
1% em etanol 
Ácido clo- Solução HCl 5 - 15 seg. Reagente universal. Permite observar detalhes muito 
rídrico 0,5 % em etanol a 85 % finos dos silicatos 
Ácido sul- Solução H, SO, 1-2 seg. Serve para diferenciar a fase intersticial, Ataca 
fúrico 1 % em etanol a 85 % de preferência os aluminatos de cálcio 
Bórax Solução 1% 5-7 min. Ataca selectivamente a alite e cal livre 
em água destilada 
Potassa Solução KOH 15 seg. Diferencia a fase intersticial 


10 % em água 
destilada 


4 — OBSERVAÇÃO MICROSCÓPICA QUALITATIVA 


O exame microscópico qualitativo dá a conhecer 
a estrutura e textura do clínquer e deixa deduzir dos 
aspectos observados quais as condições experimentais 
de cozedura e arrefecimento industriais e elucida tam- 
bém sobre o grau de moagem e homogeneização do 
cru. À miscroscopia permite portanto uma paragénese 
mineralógica, ou seja, deduzir o processo tecnológico 
de fabricação do clinquer em função da microestrutura 
do material [7, 9]. Assim, por exemplo, agrupamentos 
localizados («ninhos») de cal livre, periclase ou belite, 
denotam uma granulometria grosseira do calcário, dolo- 
mite ou quartzo, resultante em geral duma moagem 
deficiente das matérias-primas; orlas de belite sobre 
a superfície dos grãos de clíinquer podem indicar um 
enriquecimento em cinzas ácidas (siliciosas) do com- 
bustível; zonas extensas de cristais pequenos junta 
mente com cal livre são indício de temperatura insufi- 
ciente na clinquerização; cristais muito grandes e quase 
soldados entre si de alite ou belite são característicos de 
sobrecozidos; fase intersticial bem diferenciada e belite 
secundária são indicação dum arrefecimento lento do 
clíinquer; fase intersticial de cristais encastrados, sem 
belite secundária, indica, pelo contrário, um arrefeci- 
mento muito rápido (têmpera); com uma atmosfera redu- 
tora do forno aparecem nódulos de ferro metálico; etc. 

Antes do verdadeiro exame microscópico é interes- 
sante, em alguns casos, observar com uma binocular 
ou lupa as características exteriores gerais dos grãos 
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de clínquer, como forma, aspecto da superfície e cor. 
Este exame pode designar-se por macroscopia do clín- 
quer. 

Em relação à observação miscroscópica, podem ainda 
considerar-se dois aspectos fundamentais: 


— Microscopia de baixa ampliação, 100 — 300 X, 
a qual tem por objectivo obter informações sobre a poro- 
sidade e homogeneidade do clínquer e, eventualmente, 
acção das cinzas do combustível. 

— Microscopia de alta ampliação, > 300 X, me- 
diante a qual se observa a microestrutura do clinquer; 
ampliações maiores (800 - 1000 X) são em geral uti- 
lizadas apenas para o estudo de certos detalhes muito 
finos dos cristais. 


5 — DETERMINAÇÃO QUANTITATIVA DAS FASES 


A determinação quantitativa baseia-se na analo- 
gia existente entre a proporção em área ocupada pelos 
diferentes constituintes e a proporção em volume dos 
mesmos. Para que se possam considerar as percenta- 
gens das áreas como percentagens volumétricas é neces- 
sário supor um material estruturalmente homogéneo 
e com os grãos dispostos ao acaso. Assim, no caso 
de clinqueres com abundância de «ninhos» ou dis- 
tribuição dos diferentes constituintes em zonas, deve 
preparar-se uma amostra suficientemente grande. Obti- 
das as percentagens volumétricas das diferentes fases 
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multiplicam-se esses valores pelas densidades respec- 
tivas (quadro 2) e os produtos expressam-se em per- 
centagens da soma total. A porosidade, exprimir-se-á 
sempre em percentagem volumétrica, 


QUADRO 2 


Pesos específicos das fases do clinquer 


Pesos específicos 
Fases 


g/ ecm3 
C,s 3,13 
Co 8-6 3,28 
C À 3,04 
Cs AF 3,77 
CaO 3,35 
MgO 3,60 


A determinação quantitativa pode efectivar-se de 
diferentes modos: 


— Por planimetria. Esta técnica consiste em medir 
planimêtricamente as áreas ocupadas pelos diferentes 
contituintes sobre fotomicrografias. É naturalmente um 
processo demasiado moroso e dispendioso pelo elevado 
número de fotografias que exige e não tem, portanto, 
significado prático. 


— Por medição de comprimentos — método de Ro- 
siwal. O método consiste em sobrepor à imagem do clín- 
quer uma escala linear e medir os segmentos secciona- 
dos pelas diferentes fases sobre esta escala, Pode 
empregar-se uma platina integradora equipada com para- 
fusos micrométricos (tipo Leitz) [2, 20] ou uma ocular 
integradora [211]. 


— Por contagem de pontos — método de Gagalov. 
Neste método a superfície do provete, convenientemente 
atacada, é coberta com uma rede uniforme de pontos 
que acertam com as diferentes fases em casualidade 
estatística. Contam-se para cada espécie mineral os 
pontos respectivos e a operação repete-se, noutra posi- 
ção do provete, até se obter um número total de acer- 
tos, ou pontos indicadores, suficientemente grande que 
garanta uma precisão aceitável. Então, o número de 
pontos, Pi, contados para a fase Fi, dividido pelo número 
total de pontos indicadores, P, definirá o teor em volume 
dessa fase: 


Pi 
cenas: DO PORN 
P 


% vol. F 


O processo baseia-se assim em cálculos estatísti- 
cos. É especialmente útil o emprego duma ocular de 
integração com rede de pontos [9, 23, 24]. 

É este último o método actualmente utilizado no 
LNEC e que, portanto, se descreve mais em pormenor, 
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5.1 — Ocular automático de 


pontos 


integradora e contador 


Emprega-se uma ocular integradora Zeiss do tipo |, 
contendo uma rede de 25 pontos dispostos no interior 
dum círculo que representa o campo de contagem, Este 
é um pouco menor do que o campo visual da ocular, de 
modo a que o observador possa apreciar simultânea- 
mente as vizinhanças do campo de contagem (fig. 7). 


Fig. 7 —- Distribuição dos pontos na ocular integradora (Zeiss) 


À ocular é constituída por um sistema Kpl 8 X que 
permite focar, por rotação da lente ocular, a rede de 
pontos. Quando se emprega um tubo binocular deve 
montar-se no outro tubo uma ocular normal da mesma 
ampliação e igualmente corrigida. 

O princípio de aplicação da rede é simples, pro- 
cedendo-se como foi já indicado (fig. 8). 


Fig. 8 —- Princípio do método por contagem de pontos 


As contagens devem fazer-se de maneira a cobrir 
o mais equitativamente possível toda a secção do pro- 
vete, Às contagens podem simplificar-se por simples 
rotação da ocular sem deslocar o provete. No entanto 
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